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We use three-dimensional molecular dynamics simulations of magnetically interacting pancake 
vortices to study vortex matter in disordered, highly anisotropic materials such as BSCCO. We 
observe a sharp 3D-2D transition from vortex lines to decoupled pancakes as a function of relative 
interlayer coupling strength, with an accompanying large increase in the critical current remniscent 
of a second peak effect. We find that decoupled pancakes, when driven, simultaneously recouple and 
order into a crystalline-like state at high drives. We construct a dynamic phase diagram and show 
that the dynamic recoupling transition is associated with a double peak in dV/dl. 
PACS numbers; 74.60Ge, 74.60. Jg 



In highly anisotropic superconductors such as BSCCO, 
the vortex lattice is composed of individual pancake vor- 
tices 0]. These pancakes, which interact both magnet- 
ically and through Josephson coupling, align under cer- 
tain conditions into elastic lines resembling those found in 
isotropic superconductors. Three-dimensional (3D) line- 
like behavior has been observed in transformer geome- 
try measurements , muon-spin- rotation |^ and neutron 
scattering [Q. Under different conditions, however, the 
pancake vortices in each layer move independently of the 
other layers, and the system acts like a stack of indepen- 
dent thin film superconductors. Such two-dimensional 
(2D) behavior has also been seen in transformer experi- 
ments . Thus in layered superconductors two different 
effective dimensionalities of the vortex pancake lattice 
may appear, each with different characteristic properties. 

Layered superconductors exhibit a striking second 
peak in magnetization measurements corre- 
sponding to an abrupt increase in the critical current of 
the material as the applied field is increased. This sec- 
ond peak is especially sharp in BSCCO, as shown in lo- 
cal Hall probe measurements 10 1 and recent Joseph- 
son plasma frequency measurements |pd| . The lattice 
appears ordered at fields below the transition and dis- 
ordered above. There seems to be no widely accepted 
agreement on the mechanism behind this transition, al- 
though numerous scenarios have been suggested, includ- 
ing vortex entanglement JTsf , dislocation proliferations 
| p2t , dynamic effects [H^, or a 3D to 2D transition in the 
vortex pancake lattice[|[0jl|,|l|,|l|l . The effect of strong 
disorder on a possible 3D-2D transition is unclear, and 
also it is not known how the transport properties would 
be affected by a 3D-2D transition. 

In 2D systems with uncorrelated pinning, the vortex 
lattice can be dynamically reordered by an applied driv- 
ing current, passing from a glassy state at zero drive, 
through plastic flow , to a reordered state at high cur- 
rent p0|-p3|. In layered systems, when the between plane 



interactions of pancakes is weak enough that the pan- 
cakes are decoupled, the pancakes on each plane should 
reorder when a high enough driving current is applied. 
This reordering may also be accompanied by a dynam- 
ically driven recoupling transition pof , but it is unclear 
where this transition could be located in relation to the 
2D reordering transition, and how it is affected by the 
externally applied magnetic field. 

To address the issue of possible 3D-2D transitions in 
disordered anisotropic materials, we have developed a 
simulation which allows for decoupling by incorporating 
the correct magnetic interactions between vortex pan- 
cakes. We show that a sharp 3D-2D transition occurs 
when the relative strength of interlayer and intralayer 
pancake interactions is varied, and that this transition is 
associated with a sharp increase in the critical current. 
Furthermore, the system exhibits a rich array of dynamic 
3D phases when driven by a current |20|. We show that 
2D decoupled pancakes can be dynamically recoupled in a 
transition that occurs simultaneously with the dynamic 
ordering in each plane. We construct a phase diagram 
as a function of interlayer coupling and applied driving 
force, and show that a sharp, experimentally observable 
second peak in dV/dl is associated with the dynamic re- 
coupling transition. 

The magnetic interactions between pancakes in a lay- 
ered material have the same logarithmic form present in 
thin films, but are highly anisotropic We have 

performed simulations in which pancakes in all layers in- 
teract magnetically with long range interactions p5| , p6t , 
in contrast to other simulations, which treated only near- 
est layer interactions ]27| , p8[ |. Our approach compli- 
ments calculations based on Lawrence-Doniach models 
p9[ . The pancake interactions are long range both in 
and between planes, and are treated according to Ref. 
pot . Josephson coupHng is neglected as a reasonable ap- 
proximation for materials in which the anisotropy 7 is 
sufficiently large p^. 
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FIG. 1. Critical current fc (filled circles) and interlayer cor- 
relation Cz (open squares) versus interlayer coupling strength 
Sm. A transition from 2D to 3D behavior occurs between 
Sm ~ '2 and Sm ~ 2.5. Inset: Critical current /c for samples 
with L = 4, 8, 12, and 16 layers. Filled diamonds: fc for 
Sm = 0.5, in the 2D regime. Open triangles: fc for Sm = 5.0, 
in the 3D regime. 




FIG. 2. (a) Diamonds: Vx = (l/Ny) " Vx', heavy line: 
dV/dl, for a sample with Sm = 1.5 that undergoes dynamic 
recoupling. (b) Pg, the fraction of sixfold coordinated vor- 
tices, for Sm ~ 1.5. (c) Cz, the interlayer correlation, for 
Sm = 1.5. (d) Diamonds: 14; heavy line: dV/dl, for a sample 
with Sm = 4.0, which remains coupled at all drives, (e) Pe, 
the fraction of sixfold coordinated vortices, for Sm = 4.0. (f) 
Cz, the interlayer correlation, for Sm = 4.0. 



The overdamped equation of motion, at T = 0, for vor- 
tex i is given by fi = - I]j=i ViU(pi j , Zij) + f"^ + fd = 
Vi, where Ny is the number of vortices, and pi.j and 
Zjj- are the distance between pancakes i and j in cylin- 
drical coordinates. The system has periodic boundaries 
in-plane and open boundaries in the z direction. The 
magnetic energy between pancakes is [0| 



U(p,,„0) = 2deo ( (1 - + ^E, 



fexp(-z/A) In — - E2 



where Ei = dp' exp{p' / X) / p' , 

E2 = J"° dp' exp{R/\)/p', R = 22. 6A, the maximum 
radial distance, eo — ^q/{AttX)'^, X is the London pen- 
etration depth, d — 0.005A is the interlayer spacing, and 
^ is the coherence length. We vary the relative strength 
of the interlayer coupling using the prefactor . The un- 
correlated pins are modeled by parabohc traps that are 
randomly distributed in each layer. The vortex-pin inter- 



action is given by = 'Ek=i ifp/Qi'^^ ' ^k'')^ii^P " 
|ri— r[f'' |)/A), where the pin radius = 0.2A, the pinning 
force fp — 0.02/g, and /q = eo/A. The case of stronger 
pins is considered in psf . The driving current must be 
increased slowly enough for the system to equilibrate at 
each drive [^. Here fd is increased by 0.00025/q every 
35000 time steps. We have simulated a 16A x 16A system 
with a vortex density of = 0.35/A^ and a pin density 
of Up = 1.0/A^ in each of L = 8 layers. This corresponds 
to Ny = 89 vortices and Np = 256 pins per layer, with a 



total of 712 pancake vortices. We have checked for finite 
size effects on systems containing up to 42 layers and 
3738 pancakes. 

In the equilibrium state of the vortex lattice at zero 
drive, we find a sharp 3D-2D transition from vortex pan- 
cakes to vortex lines when the strength Sm of the in- 
terlayer coupling is decreased, as shown in Fig. |l|. We 
quantify the transition by measuring the spatial correla- 
tion of pancakes in neighboring planes, = 1 — ((| (r^^^ — 
ri,L-M)|/(ao/2))e(ao/2 - \{v,^l - y,^l+i)\)), where oq is 
the vortex lattice constant. A clear sharp drop in with 
decreasing Sm appears in Fig. |l| at s,„ — 2.0 for a sample 
with L = 8 layers. The 3D-2D transition is accompanied 
by a large increase in the critical current /c, as seen in 
Fig. |l|. For weak interlayer coupling, Sm < 2.0, the dif- 
ferent layers of the sample depin independently and fc is 
close to the value it would have in a 2D sample. Here, 
fc is insensitive to the number of layers in the system, 
as can be seen by comparing the data from samples with 
L =4 to 16 in the inset of Fig. |^. In contrast, coupled 
lines of pancakes at Sm > 2.0 average the random pin- 
ning over their length and become much less effectively 
pinned. Therefore fc decreases with increasing number 
of layers as seen in the inset of Fig. ^ 

When the magnetic field H increases, pancakes within 
a plane are brought closer together, but the distance 
between planes is unchanged. Thus increasing H cor- 
responds to weakening the coupling between planes by 
decreasing Sm- Therefore, our results support the sug- 
gestion that the sharp second peak observed in magneti- 
zation measurements results from a dimensional change 
in the vortex lattice from weakly pinned 3D vortex lines 
to strongly pinned 2D pancakes as the magnetic field is 
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FIG. 3. Phase diagram for varying interlayer coupling 
Sm and driving force fd- Filled circles: depinning line. X's: 
smectic transition line (in 3D phase only). Triangles: in- plane 
reordering line. Diamonds and heavy dashed line: recoupling 
line. For samples with < 2.0, a dynamic transition into 
the 3D reordered state occurs at high drives. Inset: Vor- 
tices (filled circles) and vortex trajectories (lines) at two dif- 
ferent currents, (a) and (b): Vortex trajectories in top and 
bottom layers, respectively, for a system with Sm = 1-5 at 
fd/ fo = 0.0125, in the 2D plastic flow regime. The vortex 
positions and flow are different in each layer, (c) : Vortex tra- 
jectories for a system with Sm = 1-5 at fd/ fo ~ 0.02, in the 
3D reordered flow regime. 

increased The behavior in fc that we observe in- 

dicates that a sharp change in transport properties can 
occur in a disordered system as a result of a change in 
the effective dimensionality. Furthermore, the sharpness 
of the transition found both here and in experiments sug- 
gests that the 3D-2D transition is first order. 

We next consider the question of a possible dynamic 
recoupling transition by applying a driving force. When 
the vortex lattice begins to move for /d > /c, it undergoes 
plastic tearing due to the strong pinning in our sample. 
For decoupled samples with s,„ < 2.0 each plane per- 
forms independent 2D plastic flow, as seen by the differ- 
ent vortex positions and trajectories in the top and bot- 
tom layer of the sample shown in Fig. ||(a,b). When high 
drives are applied, however, the pancakes recouple into 
lines and all planes begin to move in unison [Fig. ^(c)]. 

The reordering transition is shown in more detail in 
Fig. ||(a-c). Here, for a sample with Sm = 1-5, the re- 
coupling transition in Cz [Fig. ^c)] is sharp and occurs 
simultaneously with the in-plane reordering transition in- 
dicated by Pg [Fig. §(b)]. Furthermore, a sharp peak in 
dV/dl appears at the transition, which will be discussed 
in more detail below. In contrast, as seen in Fig. ^(d-f), 
a sample with = 4.0 that is above the static 2D-3D 
transition contains vortices that move as stiff 3D lines, 
and shows the same reordering transitions seen in previ- 
ous work on effectively two-dimensional systems [ p2[ . 

We summarize the behavior of the system in the phase 
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FIG. 4. (a) 14(1) curves for samples with different values 
of s„. From right to left, s„ =0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 
2.5, 3.0, 4.0, 5.0, 6.0, and 7.0. Samples with > 2.0 are cou- 
pled at all drives; samples with Sra < 2.0 undergo dynamic 
recoupling. (b) dV/dl curves for the same samples, from right 
to left, as listed above. Samples with Sm > 2.0 show a single 
peak in dV/dl associated with the plastic flow regime at cur- 
rents just above the 3D depinning current of fd/ fo = 0.0015. 
Samples with Sm < 2.0 show two peaks in dV/dl. (c) Filled 
squares: The value of dV/dIo, the maximum in the dV/dl 
curve, for samples with different s^. dV/dIo increases signifi- 
cantly as the transition value of Sm is approached from below. 
Open circles and dashed line: fc versus Sm. (d) dV/dl for a 
single sample with Sm = 1.0 illustrating the double peak fea- 
ture. The first broad peak (left arrow) appears just above 
the 2D depinning current of fd/ fo ~ 0.010 and is associated 
with 2D decoupled plastic fiow. The second sharp peak (right 
arrow) is associated with the dynamic recoupling transition. 

diagram of Fig. ||. The coupled vortices with Sm > 2.0 
undergo plastic flow of stiff lines, pass through a smectic 
state, and finally reorder into a crystalline-like state at 
high drives. The plastic flow and smectic regions shrink 
as the number of layers is increased, and finally disap- 
pear, so that for L = 16 and s,„ = 5 we observe elastic 
depinning directly into the 3D ordered state, with no 
plastic flow. Decoupled vortices with Sm < 2.0 depin 
into a 2D plastic flow phase in which each layer moves 
independently. The vortices switch abruptly from 2D to 
3D behavior at the recoupling transition line, so they 
directly enter the 3D ordered state. The dynamic recou- 
pling transition line and in-plane reordering transition 
line fall on top of each other in the phase diagram. 

As shown in Fig. ^(b), the dynamic recoupling and si- 
multaneous reordering are associated with a sharp peak in 
the dV/dl curve. This peak is distinct from the broader 
peak in dV/dl associated with plastic flow of the vortex 
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lattice, as indicated in Fig. ^(d). The sharp peak dis- 
appears into the background value of dV/dl and is not 
observed when the interlayer coupling becomes too weak. 
The height of the reordering peak increases rapidly as the 
static 2D-3D transition value of Sm is approached from 
below, as indicated in Fig. ^(c), and simultaneously the 
current at which the reordering peak appears shifts down- 
wards towards the location of the broad plastic peak. To 
understand the reordering peak, note that the 2D decou- 
pled lattice dcpins at the high fc associated with the 2D 
limit. When the lattice recouples, it must cross from the 
low Vx response in the 2D plastic limit up to the higher 
3D elastic Vx response curve at the recoupling transition. 
The location of the depinning transition is not affected by 
the value of Sm since it is purely a 2D effect; however, the 
recoupling transition is shifted to lower driving currents 
as Sm increases. As a result, the rising V(I) curve be- 
comes steeper as the static transition point is approached 
from below, as can be seen for Sm =1.0, 1.5, and 2.0 in 
Fig. |(a). The maximum value of dV/dl, which we call 
dV/dIo, correspondingly increases, as shown in Fig. ^(c). 
When the static 2D-3D transition is crossed at Sm = 2.0, 
the second sharp peak disappears. The behavior of this 
sharp peak in dV/dl, associated with the dynamic recou- 
pling transition, should be experimentally observable in 
transport measurements performed at fields approaching 
the second peak from above, when the vortex pancakes 
are expected to be decoupled. 

In summary, we have used a 3D molecular dynam- 
ics simulation employing the magnetic interactions of 
pancake vortices to study the dynamic phases of vortex 
matter in disordered highly anisotropic materials such 
as BSCCO. As a function of the relative interlayer cou- 
pling strength, we observe a sharp 3D-2D transition from 
vortex lines to decoupled pancakes. We find an abrupt 
large increase in the critical current as the 3D-2D line is 
crossed in a direction corresponding to increasing H, with 
decoupled pancakes being much more strongly pinned. 
At driving currents well above depinning, we find that the 
decoupled pancakes simultaneously recouple and reorder 
into a crystalline-like state at high drives. We construct 
a phase diagram as a function of interlayer coupling and 
show that the recoupling transition coincides with the 
single-layer recrystallization transition. We show that 
the recrystallization is associated with an experimen- 
tally observable double peak in dV/ dl and that the peak 
height grows rapidly as the static recoupling transition 
point is approached from below. 
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